Non-equilibrium electroweak baryogenesis at preheating after inflation 

Juan Garcia-Bellido 
Theoretical Physics, Blackett Laboratory, Imperial College, Pnnce Consort Road, London SW7 2BZ, U.K. 






s 



00 

> 

O 
(^ 

a^ 
Oh: 

CD . 



X 



Dmitri Grigoriev 
Institute for Nuclear Research of Russian Academy of Sciences, Moscow 117312, Russia 

Alexander Kusenko 
Department of Physics, University of California, Los Angeles CA 90095-1547, U.S.A. 

Mikhail Shaposhnikov 
Institute for Theoretical Physics, University of Lausanne, CH-1015 Lausanne, Switzerland 

(February 23, 1999) 

We present a novel scenario for baryogenesis in a hybrid inflation model at the electroweak scale, 
in which the Standard Model Higgs fleld triggers the end of inflation. One of the conditions for 
a successful baryogenesis, the departure from thermal equilibrium, is naturally achieved at the 
stage of preheating after inflation. The inflaton oscillations induce large occupation numbers for 
long-wavelength configurations of the Higgs and the gauge fields, which leads to a large rate of 
sphaleron transitions. We estimate this rate during the first stages of reheating and evaluate the 
amount of baryons produced due to a particular type of higher-dimensional CP violating operator. 
The universe thermalizes through fermion interactions, at a temperature below critical, Trh ^ 100 
GeV, preventing the wash-out of the produced baryon asymmetry. Numerical simulations in (1+1) 
dimensions support our theoretical analyses. 
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I. INTRODUCTION 

One of the most appealing explanations for the baryon 
asymmetry of the universe utilizes the non-perturbative 
baryon-number-violating sphaleron interactions present 
in the electroweak model at high temperatures W,0|. In 
addition to B, C and CP violating processes, a depar- 
ture from thermal equilibrium is necessary for baryogen- 
esis ^. The usual scenario invokes a strongly first-order 
phase transition to drive the primordial plasma out of 
equilibrium and set the stage for baryogenesis ||]. This 
scenario presupposes that the universe was in thermal 
equilibrium before and after the electroweak phase tran- 
sition, and far from it during the phase transition. Al- 
though there is a mounting evidence in support of the 
standard Big-Bang theory up to the nucleosynthesis tem- 
peratures, C(l) MeV, the assumption that the universe 
was in thermal equilibrium at earlier times is merely a 
result of a (plausible) theoretical extrapolation. In this 
paper we propose a picture of the early universe in which 
thermal equilibrium is maintained only up to tempera- 
tures of the order of 100 GeV. The earlier history of the 
universe is diluted by a low-scale period of inflation, after 
which the universe never reheated above the electroweak 
scale. 

We will show that the absence of the usual thermal 
phase transition at the electroweak scale does not pre- 
clude electroweak baryogenesis. In fact, according to 



recent studies of reheating after inflation, the universe 
could have undergone a period of "preheating" [Q , during 
which only certain modes are highly populated, and the 
universe remains very far from thermal equilibrium [g|. 
Such a stage creates an ideal environment in which a 
substantial baryon asymmetry could be created. The 
sphaleron transitions, known to cause a baryon number 
violation at high temperature, may also proceed in a sys- 
tem out of thermal equilibrium. In addition, the very 
non-equilibrium nature of preheating may facilitate the 
baryon number generation, as emphasized in Ref. Ig] in 
the context of GUT baryogenesis. 

It remains a challenge to construct a natural model 
with a low scale of inflation. The main problem is to 
achieve an extreme flatness of the effective potential in 
the inflaton direction (i.e. the smallness of the inflaton 
mass) without flne-tuning. Although several models have 
been proposed [R-pl, the lack of naturalness remains a 
serious problem. Perhaps recent ideas |po|-p^ related to 
large internal dimensions can provide a solution. In our 
paper we will not address the problem of naturalness, but 
will simply assume that the electroweak-scale inflation 
took place. The main question we are going to address 
is whether the electroweak baryogenesis could take place 
under these circumstances. The only qualitative feature 
of the low-energy inflation that is essential to us is that it 
produces a "cold" state in which coherent bosonic flelds 
are misplaced from their equilibrium vacuum values. An- 



other mechanism that can produce a similar state is re- 
lated to strong supercooling and spinodial decomposition 
phase transition which can occur, for example, in theories 
with radiative symmetry breaking |13| ]. 

As a toy model, we consider a hybrid model of in- 
flation 0, in which the inflaton is a SU(2)xU(l)-singlct 
and the ordinary Higgs doublet is the triggering field that 
ends inflation. Alternatively, one can view this process as 
one in which the inflaton coupling to the Higgs induces 
dynamical electroweak symmetry breaking, when the in- 
flaton slow-rolls below a certain critical value. The reso- 
nant decay of the low-energy inflaton can generate a high- 
density Higgs condensate characterized by a set of narrow 
spectral bands in momentum space with large occupation 
numbers. The system evolves towards equilibrium while 
slowly populating higher and higher momentum modes. 
The expansion of the universe at the electroweak scale is 
negligible compared to the mass scales involved, so the 
energy density is conserved, and the final reheating tem- 
perature Ti-h is determined by the energy stored initially 
in the inflaton field. We will find model parameters such 
that the final thermal state has a temperature below the 
electroweak scale, T^h < Tc ~ 100 GeV. 

Sphalerons are large extended objects sensitive mainly 
to the infrared part of the spectrum. We will conjec- 
ture that the rate of sphaleron transitions at the non- 
equilibrium stage of preheating after inflation can be es- 
timated as Fsph ^ a^Tjfjj, where Toff is some "effective" 
temperature associated with the long wavelength modes 
of the Higgs and gauge flelds that have been populated 
during preheating. 

Since Trh < Tc, the baryon-violating processes, rel- 
atively frequent in the non-thermal condensate, are 
strongly suppressed as soon as the plasma thermalizes 
via the interaction with fermions. Therefore, the baryon 
asymmetry created at the end of preheating is not washed 
out. This is in contrast to the equilibrium electroweak 
baryogenesis, where the main constraint arises from the 
tendency for the baryon density to equilibrate back to 
zero during the slow cooling following the electroweak 
phase transition. Since the energy density at the elec- 
troweak scale is so low, the universe expansion is essen- 
tially irrelevant and does not affect the prediction for the 
baryon asymmetry. 

The paper is organized as follows. In section II we 
discuss the hybrid inflationary model suitable for non- 
equilibrium electroweak baryogenesis. We estimate the 
sphaleron rates and produced baryon asymmetry in sec- 
tion HI. Our estimates, based on a number of assump- 
tions about the complicated non- linear dynamics, are in 
agreement with numerical simulations discussed in sec- 
tion IV. We summarize our conclusions in section V. 



II. AN INFLATIONARY MODEL FOR THE 

ELECTROWEAK BARYOGENESIS DURING 

PREHEATING 

Inflation is often associated with processes occurring 
at the very high energy scales, of order the Grand Unifi- 
cation Scale (~ 10^^ GeV), see Ref. |l|l. However, this 
need not be the case. Low-scale inflation models have 
been considered [0-§|. One of the side benefits of low- 
ering the inflation scale is avoiding the gravitino over- 
production constraints. There are other particle-physics 
motivations for using the TeV scale, which is associated 
with supersymmetry breaking in a class of models [[l5| . 
We will discuss a simple hybrid inflation model which sat- 
isfies the constraints from cosmic microwave background 
(CMB) anisotropics and large-scale structure, and at the 
same time contains the desired features for a successful 
reheating of the universe. We ignore completely the issue 
of radiative corrections, as discussed above. 

As described in the introduction, we want to construct 
a model with a reheating temperature which is below 
that of the electroweak scale, so that sphaleron processes 
are suppressed after reheating. Such a model necessar- 
ily has a very low rate of expansion during inflation, 
H ^ p^/"^ /Mp « 10~^ eV, which is many orders of mag- 
nitude smaller than the mass scales we will consider. This 
means that essentially all the energy density during in- 
flation is converted into radiation in less than a Hubble 
time, i.e. before the universe has had a chance to expand 
significantly. This imposes a very strict constraint on the 
energy scale during inflation^. For example, if we want 
the universe to reheat to Trh ^ 100 GeV, we need a model 
of inflation with an energy density of order p^'* ^ 200 
GeV. We will construct an example of such a model. 

Hybrid inflation M is an ingenious model of inflation, 
in which the amplitude of CMB anisotropics is not neces- 
sarily related to the GUT scale physics [|l^. The idea is 
very simple: instead of ending inflation via deviations 
from the slow-roll, it is the symmetry breaking by a 
scalar field coupled to the inflaton that triggers the end 
of inflation. The model can then satisfy the CMB con- 
straints M;f7| and allow for the electroweak-scale infla- 
tion. We will assume that the symmetry breaking fleld 
is in fact the Standard Model Higgs field, and that the 
inflaton is an additional SU(2)xU(l)-singlet scalar field. 
The model thus contains two fields, the inflaton a with 
mass m, coupled, with coupling g, to the Higgs field 
H"^ H — (/'^/2, with false vacuum energy Vb = Af'*/4A 
and the vacuum expectation value 0o — AI/\^ = v, 



^One can avoid such constraints by coupling the inflaton to 
some additional hidden-sector fields that do not contribute to 
the reheating of the observable universe. Then the potential 
energy density during infiation can be significantly larger than 
p~(200GeV)'». 
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During inflation, the inflaton is large, a ^ ac = M/g, 
and the effective mass of (p is, therefore, large and pos- 
itive. As a consequence, the Higgs field is fixed at 
(/> = and does not contribute to the metric perturba- 
tions that gave rise to the observed CMB anisotropics. 
As the inflaton field slowly rolls in the effective potential 
V{a) = Vo -I- m^tT^/2, it will generate the perturbations 
observed by COBE on large scales ||l^. Eventually, the 
inflaton reaches ct = CTc, where the Higgs has an effective 
zero mass, and at this point the quantum fluctuations of 
the Higgs fleld trigger the electroweak symmetry break- 
ing and inflation ends. The number of e-folds of inflation 
required to solve the horizon and flatness problems is 
given by 
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The fluctuations seen by COBE on the largest scales 
could have arisen in this model, N^ ~ 34 e-folds before 
the end of inflation. The observed amplitude and tilt of 
CMB temperature anisotropics 0,|l3l, ^T/T ~ 2 x 10"^ 
and n— 1 < 0.1, imposes the following constraints on the 
model parameters |18|: 
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For example, for v — 246 GeV (the electroweak symmetry 
breaking vacuum expectation value), A ~ 1, and g ~ 0.1, 
we find m ~ 2 X 10~^^ eV, and it turns out that the 
spectrum is essentially scale- invariant, n— 1 ~ 5 x 10^^"*. 
These parameters give a negligible rate of expansion dur- 
ing inflation, H ~ 7 x 10~^ eV, and a reheating temper- 
ature Trh — 70 GeV. However, the relevant masses for us 
here are those in the true vacuum, where the Higgs has 
a mass m^ = V2\v ~ 350 GeV, and the inflaton field a 
mass m — gv ~ 2^ GeV. Such a field, a singlet with re- 
spect to the standard model (SM) gauge group, could be 
detected at future colliders because of its large coupling 
to the Higgs field @. 

Some comments are in order. The consideration car- 
ried out below is qualitatively applicable also to a more 
complicated theory than the minimal SM. Let us take the 
minimal supersymmetric standard model (MSSM) with 
an additional singlet field, the inflaton cr, as an example. 
There are three SU(2) invariant couplings of the infla- 
ton to the Higgs doublets Hi and H2: ^ncr^ea/sfff -fff , 
922<^^ ^apH^ H2 ^ and gi2(7^ tapHf H^ . The Higgs mass 
matrix of the MSSM has the eigenvalues that range 
from the lightest, ~ 100 GeV, to the heaviest, roughly, 
500 GeV jl^] . In general, the inflaton-Higgs interaction is 
not diagonal in the basis that diagonalizes the Higgs mass 
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FIG. 1. The projected effective potential V{a)/Vo, for the 
inflaton field ct/ctc after the end of inflation. The dashed line 
corresponds to the m?a^ approximation around the minimum 
of the inflaton potential. Due to the shape of the potential at 
large cr, initial large-amplitude oscillations of the fleld o are 
not exactly harmonic. 



matrix in the broken-symmetry vacuum. In fact, the en- 
tire Higgs mass matrix is important in determining the 
conditions for parametric resonance. We will leave the 
analysis of multiple Higgs degrees of freedom for future 
work because it is too complicated and is not necessary 
to illustrate the main idea. 



A. Preheating in hybrid inflation 

To study the process of parametric resonance after 
the end of inflation in this model, let us recall some of 
the main features of preheating in hybrid inflation [gO[ . 
In hybrid models, after the end of inflation, the two 
fields a and (j) start to oscillate around the absolute 
minimum of the potential, cr = and (f> — v, with 
frequencies that are much greater than the rate of ex- 
pansion. Other bosonic and fermionic fields coupled to 
these may be parametrically amplified until the backre- 
action occurs and further rescattering drives the system 
to thermal equilibrium. Initially, rescattering of the long- 
wavelength modes among themselves drives them to lo- 
cal thermal equilibrium, while only a very small fraction 
of the short-wavelength modes are excited. The spec- 
tral density evolves slowly towards the higher and higher 
momenta [^,|2|. Eventually, thermalization should oc- 
cur through a process that breaks the coherence of the 
bosonic modes, e.g. through the decay of the Higgs or 
gauge fields into fermions. Such a process is very fast 
in the absence of the expansion of the universe. What 
prevented the universe from reheating immediately after 
inflation in chaotic models was the fact that the rate of 
expansion in those models was much larger than the de- 
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FIG. 2. The growth parameter /ik as a function of mo- 
menta k, in units of m, for a Higgs mass m^ = 350 GeV. The 
occupation numbers for each mode k can be obtained from 
Uk = exp(2/ifemt)/2. 



cay rate of the inflaton, and particles did not interact 
with each other until the rate of expansion dropped be- 
low the decay rate. In our case, the opposite is true: the 
rate of expansion H ~ 10~^ eV is much smaller than 
the typical gauge field decay rate into fermions, and the 
universe thermalizes quickly. Since the masses are much 
greater than the rate of expansion, many oscillations (of 
order 10^^) occur in one Hubble time |p3|. It is, there- 
fore, possible to approximate the particle production by 
that in a flat Minkowski space-time p3| . 

The evolution equation for the Fourier component of 
the Higgs field that is subject to parametric resonance is 
approximately given by 



4>k + [fc2 - Af2 + 3A(02) + g^<j^t)]cj,k = . 



(5) 



Note that this equation applies only in the case X ^ g^, 
where we have ignored the non-linear effect of the inflaton 
field a, and in particular the cross-terms g'^(f>k<^k, which 
do not contribute significantly before backreaction [ pO[ . 
We will only use this equation for qualitative arguments, 
since our quantitative results will be fully non-linear and 
non-perturbative, based on numerical simulations, see 
Section IV. As the inflaton oscillates around a = with 
amplitude "E = ac = M/g in the effective potential of 
Fig. 1, its coupling to the Higgs will induce the paramet- 
ric resonance with a q parameter [M , characterizing the 
strength of the resonance, and given by 
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Since we can neglect the rate of expansion, the ampli- 
tude of oscillations S does not decrease, and the reso- 
nance is extremely long-lived. For generic values of the 
couphngs, g^ ~ 10~^ — 10"'^, it is, in fact, a broad res- 
onance, q :s> 1. Higgs particle production occurs at the 
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FIG. 3. The evolution of the Higgs spectrum UkUDk, in 
units of « = 246 GeV, from time to 10* v~^ , as a function of 
momentum, k/m. The initial spectrum is set by preheating, 
and contains a set of narrow bands (solid line). The subse- 
quent evolution of the system leads to a redistribution of en- 
ergy between different modes. Note how rapidly a "thermal" 
equidistribution is reached for the long-wavelength modes. 
However, the whole Higgs spectrum approaches thermaliza- 
tion already in the middle of resonance (see Fig. W below). 

instants when cr(t) = 0, and continues until backreac- 
tion becomes important, cither for the inflaton oscilla- 
tions {{(jP') ~ m^/g'^) [|l]j2|,|l or for the effective Higgs 
mass (3A(02) ~ M^ -~^+ 4m^y/q) [|. Which of the 
two effects back-reacts first depends on the coupling g. 
Here k^ — \/2m q^l"^ is the typical momentum of the res- 
onance band. For g > 0.08 backreaction on the inflaton 
mass occurs before the A-term in (||) is relevant. We have 
chosen g = 0.1 for definiteness, and computed the power 
spectrum of the Higgs field. For a smaller coupling g, the 
resonance spectrum would be different, but the qualita- 
tive behavior would be similar. In fact, it does not matter 
how many bands the parametric resonance populates be- 
cause after rescattering all those bands smooth out and 
reach "thcrmalization" over a finite region in momentum 
space [py,p2| . In Fig. we show the growth parameter [ik 
as a function of fc. The typical momentum contributing 
to the power spectrum, fc^l^^p, is 



k^k^/2 = rn q^'^/V2 w 2m , 



(7) 



where the growth factor has a large value, /.tmax — 0.9. 
This unusually large number is due to the fact that what 
drives the Higgs production in this model is not the usual 
parametric resonance from oscillations around the min- 
imum of the potential pfl], but the spinodal instability 
responsible for the breaking of the electroweak symme- 
try. In the language of Mathieu equations p4| , this cor- 
responds to a large and negative A = (fc^ — M'^)/Am? 
parameter, which induces large growth factors /i. On the 
other hand, the occupation number of a given mode is 
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FIG. 4. The time evolution of the inflaton energy, the 
Higgs energy and the total energy. Note that the energy is 
measured in units of v and time in units of v~^ , see Ref. fesi. 



determined from ^k as n^ — i exp(2/^fcTOi). This means 
that, within a few oscillations, the Higgs field reaches a 
huge occupation number over a range of narrow bands 
in momentum space. Therefore, the Higgs fluctuations 
grow exponentially with time. 
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with /x = /ieff — 0.8. At backreaction, the Higgs ex- 
pectation value is just of order its vacuum expectation 
value (VEV), {(fr^) < m? /g^ ~ i;^, but continues to grow 
slightly during rescattering |^,^ . With our set of pa- 
rameters, this happens at times t ~ C(l) GeV~^. 

In Section IV of this paper we follow a numerical ap- 
proach in (1-1-1) dimensions and computed the initial 
state from parametric resonance and subsequent stages 
like rescattering and backreaction directly through the 
real time evolution of the classical equations of motion for 
the bosonic modes with arbitrary fc, with all the couplings 
between fields properly taken into account. This way, we 
have automatically included rescattering and thermaliza- 
tion in the evolution. 



B. Higgs coupling to W bosons 

Soon after production, Higgs particles decay predom- 
inantly into W bosons with a branching ratio of order 
one, for m^ = 350 GeV, and a decay rate F ~ 20 GeV. 
One may ask whether the Higgs oscillations may induce 
a resonant production of gauge bosons. It turns out that 
the corresponding resonance is very narrow and insuf- 
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gl^^^/Am^ ~ 0.3 GeV < F, where 



g^ = A-Ka^ is the SU(2) gauge coupling and $ ~ w/10 
is the amplitude of the Higgs oscillations during the first 
resonance stage, see Fig. 9) for the coherent decay of the 
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FIG. 5. The time evolution of the effective temperature, 
in units of v. We have averaged the Higgs power spectrum 
over different low-momentum regions, and we obtain several 
effective temperatures that show different time behavior. 



Higgs into gauge bosons. It is therefore appropriate to 
use perturbation theory to calculate the Higgs decay into 
the W bosons. 

Since the rate of growth of the energy density of the 
Higgs field. 
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is larger than its decay rate into W bosons, i.e. 2^m ~ 
2F ^ 40 GeV, we do not expect a significant depletion of 
the energy density of the Higgs field during preheating, 
while the energy density of the gauge bosons grows expo- 
nentially at the same rate, p„ oc exp(2/i7TT,i). Therefore, 
soon after rescattering, most of the energy density is in 
the form of Higgs and gauge fields with essentially zero 
momentum. It is these long-wavelength gauge configu- 
rations that will play an important role in inducing the 
sphaleron transitions, and the subsequent baryon pro- 
duction. 

One of the most fascinating properties of rescatter- 
ing after preheating is that the long-wavelength part of 
the spectrum soon reaches some kind of local equilib- 
rium ||2l| , p2| , while the energy density is drained, through 
rescattering and excitations, into the higher frequency 
modes. Therefore, initially the low energy modes reach 
"thermalization" at a higher effective "temperature" [pj , 
while the high energy modes remain unpopulated, and 
the systenr is still far from true thermal equilibrium: 
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It is possible to estimate the effective "temperature" T^s 
from the conservation of energy during preheating. The 
energy per (long wavelength) mode is UkUJk ~ T^s, or 
effectively equipartitioned. Since only the modes in the 
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FIG. 6. The variance of the Chern-Simons number, i.e. 
(iV^g)— {Nf^g(t))'^, as a function of time. The solid line is a 
result of smoothing out the measured values (dots) . 



range < /c ^ k^ax — 4fc* ~ 5mq^'^ are populated, 
we can integrate the energy density in Higgs and gauge 
fields, 5j3 = 1 + 3 X 3 = 10, to give, in (3+l)-dimensions, 
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which gives Tcff ~ 350 GeV. We note that the effective 
temperature depends on the value of the coupling g as 
T^g oc g~^, as expected [Q. The temperature Toff is 
higher than the final reheating temperature Tj-h, which 
is easy to understand, since preheating is a very efficient 
mechanism for populating just a few modes, into which 
a large fraction of the original inflaton energy density is 
put. This means that a few modes carry a large amount 
of energy as they come into partial equilibrium among 
themselves, and thus the effective "temperature" is high. 
However, when the system reaches a full thermal equilib- 
rium, the same energy is distributed between all modes, 
which corresponds to a much lower temperature. In our 
example, thermalization of long-wavelength modes hap- 
pens at a time scale t ^ T^^ ^ C'(l) GeV~^, where 
r„ ~ 2 GeV is the width of the vector boson. 

This is the main reason why the out of equilibrium 
mechanism of preheating is so efficient in producing 
sphalcron transitions, since the rate of these transitions is 
greatly enhanced by the higher effective temperature. An 
alternative way of seeing this is by analogy with a diffus- 
ing plasma. The rescattering of Higgs and W bosons after 
preheating produces a diffusion which enhances over-the- 
barrier sphaleron transitions. It is fortunate that the de- 
scription of this diffusion mechanism can be done with 
the use of an effective temperature, for which the rate of 
sphaleron transitions can be estimated analytically, by 
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FIG. 7. The sphaleron transition rate rsph(i) is propor- 
tional to the time derivative of the Chern-Simons variance 
(solid line on Fig. |6|). The left scale measures the transitions 
per unit time per whole sample volume. The right scale mea- 
sures the sphaleron transitions per unit time per unit volume. 
Note that because of high final equilibrium temperature the 
transitions aren't vanishing after the resonance. 



ignoring the higher momentum modes and the integra- 
tion over hard thermal loops pa . 



III. BARYON ASYMMETRY OF THE UNIVERSE 

It is well known that sphaleron transitions are mainly 
sensitive to the long- wavelength modes in a plasma. This 
is because the sphaleron size, (a^roff)^^, is much larger 
than the typical Compton wavelengths of particles in the 
plasma, {2k.^,)~^ ~ {5m)~^. A simple argument then sug- 
gests that the rate of sphaleron transitions per unit time 
per unit volume should be of the order of the fourth power 
of the magnetic screening length in the plasma ||2^,^ .[| 

We, therefore, conjecture that the sphaleron transition 
rate during rescattering after preheating, Fsph, can be 
approximated by that of a system in thermal equilibrium 
at some temperature Tcs defined in the previous section: 



■^In the symmetric phase of the electroweak theory and in 
thermal equilibrium at a temperature T, the typical momen- 
tum scale of sphaleron processes is a„T (a„ ~ 1/29 is the 
weak gauge coupling) which is much smaller than the aver- 
age momentum of the particles in the plasma, fc ~ T. It 
was argued in Refs. [ p9|j2(| that the higher momentum modes 
with typical scale greater than g^^T should slow down the 
sphaleron processes by an extra factor a^ log(l/a.^^). Dur- 
ing the first stages of reheating those high frequency modes 
are not populated and therefore should not be considered in 
our estimate. 
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FIG. 8. The expectation value of (a^) in units of v^, as a 
function of time. 
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In the Standard Model, baryon and lepton numbers are 
not conserved because of the non-perturbative processes 
that involve the chiral anomaly: 



d,r, = d,f^ = 



327r2 ^^"^ ■ 



(13) 



Furthermore, the sphaleron configurations connect vacua 
with different Chcrn-Simons numbers, N^^, and induce 
the corresponding changes in the baryon and lepton num- 
ber, AS = AL = SAiVcs . 

A baryon asymmetry can be generated by sphaleron 
transitions in the presence of C and CP violation. There 
are several possible sources of CP violation at the elec- 
troweak scale. The only one confirmed experimentally 
is due to Cabibbo-Kobayashi-Maskawa mixing of quarks 
that introduces some violation of CP, but it is probably 
too small to cause a sufficient baryon asymmetry. Vari- 
ous extensions of the Standard Model contain additional 
scalars (e.g. extra Higgs doublets, squarks, sleptons, etc.) 
with irremovable complex phases that lead to C and CP 
violation. 

We are going to model the effects of CP violation in the 
effective field theory approach. Namely, we assume that, 
after all degrees of freedom except the gauge fields, the 
Higgs, and the inflaton are integrated out, the effective 
Lagrangian contains some non-renormalizable operators 
that break CP. The lowest, dimension-six operator of this 
sort is M 



O 



now 



327r- 



F.^F^-^ , 



(14) 



The dimensionless parameter 5^p is an effective mea- 
sure of CP violation, and Mnow characterizes the scale 
at which the new physics, responsible for this effective 
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FIG. 9. The expectation value of {4>''4>), in units of v^ , as 
a function of time. The variations in {cj>''(j)) give rise to a CP 
nonconserving chemical potential /^cff, which is the source of 
the baryon asymmetry. 



operator, is important. Of course, other types of CP vio- 
lating operators are possible although, qualitatively, they 
lead to the same picture. 

Note that the operator (|4|) is CP-odd but does not vi- 
olate C. Thus, in a pure bosonic theory non-equilibrium 
evolution can only produce parity-odd or CP-odd config- 
urations, but no C asymmetry. For example, the Chcrn- 
Simons number can be produced, as it is C-even but 
P and CP-odd. C violation, necessary for baryogene- 
sis, conies from ordinary gauge-fermion electroweak in- 
teractions that violate C and parity, but conserve CP. 
This manifests itself in the anomaly equation that re- 
lates baryon number (C-odd but P-even operator) to the 
Chern-Simons number (C= -1-1, P= —1). In other words, 
C violation in the bosonic sector of the theory is not re- 
quired as long as it appears in the fermionic sector, via 
the electroweak interactions. 

If the scalar field is time-dependent, the vacua with dif- 
ferent Chern-Simons numbers are not degenerate. This 
can be described quantitatively in terms of an effective 
chemical potential, /ioff , which introduces a bias between 
baryons and antibaryons, 



Moff 



Af 2 ,„ dt 



(15) 



This equation follows from Eq. (O) by integration by 
parts. Although the system is very far from thermal equi- 
librium, we will assume that the evolution of the baryon 
number n^ can be described by a Boltzmann-like equa- 
tion, where only the long-wavelength modes contribute. 
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(16) 



where F^ = (39/2)Fsph/T3g. The temperature T^s de- 
creases with time because of rescattering: the energy 
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FIG. 10. The Chern-Simons number A^'^g (i) averaged over 
100 runs, for CP violating parameter k — —1. Note that it 
settles at a large and negative A'^^g , corresponding to large 
baryon production. 



stored in the low-frequency modes is transferred to the 
high-momentum modes. 

The rate F^ , even at high effective temperatures, is 
smaller than other typical scales in the problem. Indeed, 
for Tcft ~ 400 GeV, F^ ~ 0.01 GeV, which is smaU com- 
pared to the rate of the resonant growth of the Higgs 
condensate (2fim ~ 40 GeV). It is also much smaller 
than the decay rate of the Higgs into W's and the rate 
of W decays into light fermions. Therefore, the last term 
in Eq. (|lq ) never dominates during preheating and the 
final baryon asymmetry can be obtained by integrating 



dtTsphit) 



^J■ctf{t) 
Toff(i) 



F, ^- 



sph 



^eff Ml 



(17) 



where all quantities are taken at the time of thermaliza- 
tion. This corresponds to a baryon asymmetry 



45Qt^cP 



2^2 ff* Ml 



m'- 



(18) 



where g* — g^ + (7/8)5p ~ 10^ is the number of effec- 
tive degrees of freedom that contribute to the entropy 
density s at the electroweak scale. Taking (0^) ~ v^ — 
(246 GeV) 2, the scale of new physics M„ew ^ 1 TeV, the 
coupling a^ ~ 1/29, the temperatures Toft — 350 GeV 
and Trh ^ 70 GeV, we find 



FIG. 11. The variance of the topological Chern-Simons 
number (TV^g)— {N^g{t))^ for a set of runs with total en erg y 



reduced by a factor of 4, as compared to runs in Figs. M-h 

IV. NUMERICAL SIMULATIONS IN (1+1) 
DIMENSIONS 

The theoretical analysis presented above was based on 
the conjecture that the sphaleron transition rate can be 
described in terms of the effective "temperature" Toff as 
in equation (jl^). This assumption is based on the rea- 
soning given above and seems quite plausible. We have 
also verified the validity of such description in the (1-1-1)- 
dimensional numerical simulations. 

For simplicity, we consider an Abelian Higgs model in 
(1-1-1) dimensions, which was successfully used before for 
studying physics relevant to baryogenesis [p5| , pi|j3^ ] . The 
Lagrangian comprises two scalar fields and a U(l) gauge 
field: 
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where T)^ 
and 



^ da - ieA, 
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J, with e the U(l) gauge cou- 
pling, and e^^ is the totally antisymmetric tensor in 
(1-1-1) dimensions. Here CP violation is induced via the 
K4>*(j)e^i,F^^ term, which violates both C and CP. Fur- 
thermore, in (1-1-1) dimensions, the analogue of the chiral 
anomaly is the anomalous non-conservation of the gauge 
invariant fermionic current, jp = tp^^ip. 



^- ~ 3 X 10" 

s 



AP 



(^)'^2x 10-^5,,, (19) 



consistent with observations for J^p ~ 10"'^, which is a 
reasonable value from the point of view of particle physics 
beyond the Standard Model. Therefore, baryogenesis at 
preheating can be very efficient in the presence of CP 
violation that comes from new physics at Afnow ^ 1 TeV. 



^-^■^ = "4^ 



p^. 



which serves as a source of B violation. 

The corresponding equations of motion are: 
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FIG. 12. The time evolution of effective temperature Tea 
in units of v, for the reduced-energy runs. Note the smoother 
rise and decline of the effective temperature with time. 



j7 — i{(l)*d'^(j) — (j)d^^<f)*) is the charged current of 
the Higgs field. 

The numerical simulations track the real-time evolu- 
tion of the classical field configurations. The initial condi- 
tions are set by preheating as a set of narrow bands in the 
Higgs power spectrum (Fig. |^). The real-time evolution 
leads to a gradual redistribution of energy between differ- 
ent modes, including the inflaton field a itself. Note that 
the Higgs field takes a very long time to reach its VEV, 
see Fig ^ due to the presence of the CP violating term in 
Eq. (p2) , which leads to the production of baryon number 
before equilibrium. In fact, the full thermalization of the 
system should take a very long time |2^,^, while some 
other processes (e.g. interaction with fermions from the 
decay of gauge fields and/or Higgs) will induce thermal- 
ization via decoherence long before that. Thus, although 
it is technically possible to reach complete thermalization 
of the whole system including the infiaton, our simula- 
tions are necessarily limited to the vector and Higgs decay 
time scale, of the order of 50 inflaton oscillations, as in 
Fig. |. 



A. The sphaleron transition rate and the effective 
temperature 

As expected, the resonant inflaton decay quickly leads 
to a population of the long-wavelength modes of the 
Higgs field, see Fig. 0. This happens after only a few oscil- 
lations of the inflaton. At this point the long-wavelength 
modes contain a very large fraction of the total energy, 
and that leads to a noticeable increase in Teff at the be- 
ginning of the resonance, see Fig. pi 

The sphaleron transitions immediately set in. We mon- 
itor them both by calculating the Chern-Simons number, 
N^^ = J^idx^ in the temporal gauge ^o = 0, and also 
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FIG. 13. The evolution of the sphaleron rate with time for 
the reduced-energy runs. The decrease of the final effective 
temperature in Fig. [13 results in the freezing out the equilib- 
rium transitions soon after the resonance. 



by keeping track of the U(l) winding number of Higgs 
field. (Actually, no statistically significant difference be- 
tween these two quantities was observed) . To get a quan- 
titative estimate of the transition rate we measure the 
variance of N^^, i.e. 5'^ = (N^^) — {N^^{t))'^, over an 
ensemble of 100 independent runs starting from different 
field configurations that have the same energy spectrum 
as shown in Fig. 0. The preparation of the initial configu- 
ration and other peculiarities of the numerical procedure 
will be discussed in detail in a future publication. 

The variance of N^^ is shown in Fig. |6[ Its time deriva- 
tive dS'^/dt = Volume x Tgph is plotted in Fig. |^. Note 
that this relation comes naturally from the diffusion of 
the Chern-Simons number, {N^^) ~ Volume x Fsphi, 
which follows a typical Brownian motion ||28|l . Note that 
for initial parameters chosen as in Figs. 



-IIOI, the rate 

actually increases during the thermalization of the Higgs 
field. However, for our purposes, it is important that we 
get a substantial amount of sphaleron transitions right 
after the beginning of the resonance. One could slow 
down the after-resonance transitions by decreasing the 
total energy of the system by a factor of 4 (see Figs. ^l|- 
p^ below). However, this decreases the net generated 
asymmetry considerably, due to a decrease of the effec- 
tive temperature Tes, see Fig. ^ and the subsequent 
decrease in the sphaleron rate just after the resonance, 
see Fig. ^ These two sets of figures helps us gain intu- 
ition about the process of baryogenesis during preheating 
in (1+1) dimensions. 



B. The generation of the baryon asymmetry 



As is clear from Eqs. (|20|) and (^2|), the chemical po- 
tential /ioff oc ~Kdo{(j)*(f>) is non-zero only during the 



100 




3000 



FIG. 14. The time evolution ol (a^), in units ol v^ , for 
the reduced-energy runs of Figs. 11-13. 



resonance. The energy transfer from the inflaton to the 
Higgs field results in a steady shift of (</)*</)) expectation 
value, see Fig. g. This shift in VEV acts as a chemical 
potential and drives the baryon asymmetry. 

The baryon asymmetry generated by the non-equili- 
brium sphaleron transitions in the presence of a CP vi- 
olating chemical potential /icg, see Eq. (Eq), is observed 
as a non-zero value of (N^^) averaged over a computer- 
generated ensemble. As shown in Figs. ^ and |l5|, (-/V^s) 
steadily increases and eventually freezes when the expec- 
tation value (0*0) approaches a constant value and the 
chemical potential ( |l5| ) vanishes. In the early universe, 
the drift of (N^^) is eventually interrupted by the decay 
of the vector and Higgs fields into fermionsjj This leads 
to thermalization and, as long as the reheat temperature 
is sufhciently low, there is no further wash-out of the 
baryon asymmetry. We note in passing that, in our nu- 
merical simulations, the Chern-Simons number attained 
at the end (i.e. the final baryon number) is approximately 
linearly dependent on the CP violating parameter k, and, 
therefore, our estimate can be extrapolated to very small 
values of k. 



V. CONCLUSION 

There is no empirical evidence that a thermal elec- 
troweak phase transition took place in the early universe. 
However, since the only well-established source of baryon 



^The Higgs and vector decays into fermions are not included 
in our (l-l-l)-dimensional simulations. There has been recent 



progress |B4| in introducing fermions in (1-1-1) lattice simula- 
tions, but we will leave for future work such developments. 
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FIG. 15. The time evolution of {(f)* (j)), in units of u^, 
(7V^g(i)), for the reduced- energy runs of Figs. 11-13. 



and 



number non-conservation is the gauge sector of the Stan- 
dard Model, one could argue that electroweak baryogen- 
esis [1] is the only explanation for the baryon asymmetry 
of the universe that does not invoke any unknown B- 
violating new physics. This reasoning would favor the 
usual electroweak phase transition, followed by the elec- 
troweak baryogenesis, on aesthetical grounds. 

In this paper we have presented an appealing alter- 
native. We have shown that a new kind of electroweak 
baryogenesis, which still uses only the known sources of 
baryon number violation, is possible even if the reheat 
temperature after inflation was too low for a thermal 
restoration of the SU(2)xU(l) gauge symmetry. More- 
over, the departure from thermal equilibrium, necessary 
for generating a non-zero baryon number density, is nat- 
urally achieved at preheating after an electroweak-scale 
inflation. Sphaleron transitions take place during pre- 
heating, before the thermalization of the plasma. The 
baryon asymmetry can be generated through sphalerons 
in a manner similar to the usual electroweak baryogene- 
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sis 0. When the universe reaches thermal equihbrium, 
the temperature can be small enough to suppress further 
baryon-violating processes, so that the baryon asymme- 
try is not washed out. 
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NOTE ADDED 

After our paper was finished, we learned about a re- 
cent paper |g^ that also discussed baryogenesis after an 
electroweak-scale inflation. 
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